Interface solitons in thermal nonlinear media 
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We demonstrate the existence of fundamental and dipole interface solitons in one-dimensional 
thermal nonlinear media with a step in linear refractive index. Fundamental interface solitons are 
found to be always stable and the stability of dipole interface solitons depends on the difference in 
linear refractive index. The mass center of interface solitons always locates in the side with higher 
refractive index. Two intensity peaks of dipole interface solitons are unequal except some specific 
conditions, which is different from their counterparts in uniform thermal nonlinear media. 
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Nonlocal solitons have been found in many physical 
systems, such asphotorefractive crystals nematic liq- 
uid crystals @, 0], lead glasses SH, atomic vapors Q, 
and Bose- Einstein condensates [3, Q, etc. In nonlocal 
nonlinear media, various types of solitons have been stud- 
ied theoretically and experiment ally , such as vortex soli- 
tons [9], multipole solitons [HI UK, L a guerre and Her- 
mite soliton clusters 12], Ince-Gaussian solitons [l3j]. For 
nonlocal solitons, there are many interesting properties, 
for instance, large phase shift [lj], attraction between 
two dark solitons [lg, [l6[ , self- induced fractional Fourier 
transform [l7| . etc. 

Surface waves are localized waves residing at the inter- 
face between two media with different optical properties. 
They have been used to study the surface properties of 
media in physics, chemistry, and biology. In the presence 
of nonlinearity, some kinds of surface solitons have been 
found theoretically and experimentally. Recently, nonlo- 
cal surface solitons have been studied both theoretically 
and experimentally [Tcl - I^ . It is proved that the nonlo- 
cality can support various types of surface solitons, such 
as m ultip ole surface solitons [l9l - [2lj | . vortex surface soli- 
tons |l9| |. and incoherent surface solitons (22|. Surface 
dipoles, vortices, and bound states of vortex solitons are 
found to be stable at two dimensional interfaces [l9|. In 
one dimension case, multipole surface solitons are stable 
when the number of poles is less than three, whereas the 
higher-order solitons can be stable in layered thermal me- 
dia [20j]. When the interface is formed by two nonlocal 
nonlinear media, both fundamental and dipole interface 
solitons are found to be stable in the presence of the op- 
tical lattice, but dipole interface solitons do not exist in 
uniform latticeless media [21( . Here we study the surface 
waves at the interface formed by two thermal nonlinear 
media, in which the boundary force effect [23|, [24[ can 
support the existence of dipole interface solitons. 

In this paper, we demonstrate the existence of interface 
solitons in the thermal nonlinear media which own two 
different linear refractive indices. Stability analysis of 
interface solitons arc carried out. It is found that funda- 



mental interface solitons are always stable and the stabil- 
ity of dipole interface solitons depends on the refractive 
index difference between two media. Two intensity peaks 
of dipole interface solitons are unequal except some spe- 
cific conditions, which is different from their counterparts 
in uniform thermal media. 

We consider a (l+l)dimensional thermal sample occu- 
pying the region —L<x<L. The sample is separated 
into two parts at the center (x = 0). Two boundaries 
(x = ±L) and the interface are thermally conductive. 
All parameters of the two parts, such as thermal con- 
ductivity coefficient, absorption coefficient and thermal 
coefficient, are the same except the linear refractive in- 
dex. The propagation of a transverse-electric (TE) polar- 
ized laser beam is governed by the dimensionless nonlocal 
nonlinear Schrodinger equation 
(i) in the left, i.e. —L < x < 
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(ii) in the right, i.e. < x < L 



.dq 

dz 



1 d 2 g 



+ nq - n d q = 0, 



d 2 n 
dx 2 



(1) 



(2) 



Corresponding author's email address: huwei@scnu.edu.cn 



where x and z stand for the transverse and longitudinal 
coordinates scaled to the beam width and the diffraction 
length, q is the complex amplitude of optical field, n is 
the nonlinear refractive index change, and n ( i > is pro- 
portional to the linear refractive index difference between 
two media. 

The boundary conditions can be described as q(±L) = 
when the laser beam is narrow and far from two 
boundaries. The continuity conditions at the interface 
(x = 0) depend on the polarization. For the TE polarized 
wave, the continuity conditions for the transverse field 
are q{-0) = q{+0) and dq/dx\ x =- = dq/dx\ x = +0 [Hj]. 
For the nonlinear refractive index, we have the boundary 
conditions n(±L) — because two boundaries are ther- 
mally stabilized by means of external heat sinks [TH, [2(| ■ 
Since the interface is also thermally conductive, we have 
the continuity relation n(— 0) = n(+0). From the conti- 
nuity of q and the nonlinear equations in Eqs. ([1} and 
@, the derivative of nonlinear refractive index is also 
continuous. 
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We search for the soliton solutions for Eqs. ([I) and 
@ in the form q(x, z) — w(x) exp(ibz), where w(x) is a 
real function, b is the propagation constant. An iterative 
method is used to get numerical solutions for different rid 
and b. The results for fundamental and dipole interface 
solitons are shown in Figs. [Hand [2 respectively. 

In order to elucidate the stability of interface solitons, 
we search for the perturbed solutions of Eqs. (HJ and @ 
in the form q — (w + u + iv) exp(ifez), where u(x, z) and 
v(x, z) are the real and the imaginary parts of small per- 
turbations. The perturbation can grow with a complex 
rate a upon propagation. Substituting the perturbed 
soliton solution into Eqs. ([T]) and ([2]), and using the 
stationary soliton solution w(x), one can get the linear 
eigenvalue problem around stationary solution w(x), 
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(-L<x< 0), (3) 
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(0 < x < L), 



bu 



ritju + wAn, 
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where An = —2 J_ L G(x, x')w{x')u{x')dx' is the refrac- 
tive index perturbation, the response function G(x, x') — 
(x + L)(x' - L)/(2L) for x < x' and G(x,x') = (x' + 
L)(x — L)/(2L) for x > x', a r (real part of a) represents 
the instability growth rate. 

The eigenvalue problem of Eqs. ^ and (HJ) has been 
solved numerically. We find that the fundamental inter- 
face solitons are always stable in their whole domain for 
both small and large rid- For comparison, we know that 
an optical beam can form a stable fundamental soliton 
in bulk thermal nonlocal media [2M [26| , and fundamen- 
tal surface solitons are also stable in thermal nonlocal 
media [2pj . 

First, we discuss fundamental interface solitons as 
shown in Fig. [T] In our model, when n<j approaches 
zero, the interface soliton reduces to a bulk soliton. As rid 
increases, the stable interface soliton shifts itself to the 
higher index part. The soliton mass center, defined as 



= 11 



x\q\ 



2 dx/ \q\ 2 dx, will always locate in the 



left part where the linear refractive index is higher. For 
example, for rid = 0-05 as shown in Fig. [TJa), x g = —3.8 
and most energy of the soliton resides in the left part, 
while a little energy resides in the right. \x g \ increases 
monotonically as rid increases, as shown in Fig. [ljc). 
When rid > 1, almost all energy of the soliton resides 
in the left part [Fig. Hfb)]. The intensity profiles for 
larger rid are similar to those of surface solitons [H, [2(| , 
whereas the nonlinear index changes are different. For 
surface solitonSjthere is a step in the nonlinear index at 
the interface [3, H3] ; but the nonlinear index distribu- 
tion is continuous at the interface in our model. 



For much larger rid, x g and wq approach to certain val- 
ues, i.e. x g — > —6 and wo — > 4.7 for b = 0.5 and L = 30 as 
shown in Figs.QJc) andQJd). Because the boundary force 
effect [23, |24j always pushes the soliton to the center of 
the sample, the interface soliton can not shift away from 
the interface. It is obvious that the asymptotic value of 
x g is proportional to the asymptotic beam width when 
rid approaches to infinity. The asymptotic value of \x g \ 
decreases as the propagation constant b increases [same 
as the beam width shown in Fig. [He)], but increases 
as the sample size L increases [Fig. He)]. As shown in 
Fig- UXf), the energy flow, defined as U — J™ \q\ 2 dx, 
is a linearly growing function of b for different rid- It is 
worthy to note that the change of rid does not influence 
the energy flow of the fundamental interface solitons be- 
cause the peak intensity increases while the beam width 
decreases. 
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FIG. 1. (Color online) Profiles of fundamental interface soli- 
tons at (a) rid = 0.05, b — 0.5 and (b) rid — 5, b — 0.5. Solid 
red lines stand for optical intensity while dashed green lines 
stand for nonlinear refractive index. Dash-dotted line indi- 
cates the interface and L — 30. (c) Soliton mass center versus 
rid at 6 = 0.5 for different sample sizes, (d) Beam width ver- 
sus rid at b = 0.5 for different sample sizes, (e) Beam width 
versus propagation constant at different rid for L — 30. (f) 
Energy flows U versus propagation constant at different rid 
for L = 30. 

The results of dipole interface solitons are shown in 
Figs. [5] and [3] The change of the two intensity peaks of 
interface dipoles is interesting. When rid approaches to 
zero, similar to a bulk dipole soliton, the two intensity 
peaks are equal. As rid increases, the two peaks increase 
at different rates and the right peak becomes higher than 
the left one [Figs. [U[a) and^b)], whereas the mass cen- 
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FIG. 2. (Color online) Profiles of dipole interface solitons 
at (a) ra d = 0.02, (c) n d = 0.029, (e) n d = 0.05, and (g) 
rid = 1. Solid red lines stand for optical intensity while dashed 
green lines stand for nonlinear refractive index, (b), (d), (f) 
and (h) denote the propagations of dipole interface solitons 
corresponding to (a), (c), (e) and (f), respectively. Note the 
dipole interface soliton at rid = 1 in (g) and (h) is unstable. 
In all cases b = 0.5, L — 30. 



ter moves toward the left [Fig. |3[a)]. Then the left peak 
increases more quickly, and the two peaks become equal 
again when rid = 0.029 [Figs. [He) andHJd)]. Increasing 
rid sequentially, the left peak becomes higher than the 
right one and a significant part of the right peak resides 
in the left part [Figs. HJe) and OUT)]. When n d > 0.5, 
such as rid — 1 shown in Figs. [2jg) and[2jh), almost all 
energy moves into the left part and the intensity profiles 
are similar to those of dipole surface solitons, however, 
their nonlinear refractive index distributions are differ- 
ent [2fj| . Similar to fundamental interface solitons, the 
soliton center moves slowly toward the left and the beam 
width decreases monotonically as rid increases [Figs. E2a) 
and[3fb)]. In Fig. (3Jc) andEJd), the beam width de- 
creases monotonically with increasing b, whicrcas the en- 
ergy flows of dipole interface solitons increase monoton- 
ically with increasing b. The energy flows for different 
rid are almost equal for the same reason mentioned in 
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FIG. 3. (Color online) (a) Dipole mass center versus rid at 
b = 0.5 for different sample sizes, (b) Beam width versus 
rid at b — 0.5 for different sample sizes, (c) Beam width 
versus propagation constant at different n d for L = 30. (d) 
Energy flows U versus propagation constant at different rid 
for L = 30. 
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FIG. 4. (Color online) Real part of the perturbation growth 
rate for dipole interface solitons at (a) rid = 0.2, 0.1, 0.05, and 
(b) rid = 1, 5. In those cases L — 30. (c) Critical value of 
propagation constant versus rid for different sample sizes. 

An important result of this paper is that the stabil- 
ity of dipole interface solitons depends on the index dif- 
ference rid- For comparison, dipole surface solitons in 
thermal nonlocal media are stable [2(J, but dipole inter- 
face solitons in nonlocal nonlinear media with a finite 
range of nonlocality are unstable [2l| . Figure |4] presents 
the results of the stability analysis for dipole interface 
solitons. For a given rid, there exists a critical propa- 
gation constant & c , and dipoles are stable when b > b c . 
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For iid — 0.2 and 0.1, the stable regions are b > 1.8 
and b > 0.5, respectively [Fig. Eta)]. The relation be- 
tween critical propagation constant b c and rid is shown 
in Fig. Stc) for different sample sizes L. If rid is very 
small, the dipole interface solitons are stable almost in 
their whole domain, which consists with the dipole soli- 
tons in bulk nonlocal media [25|, [26| . 

In Ref. [ll|, lattices are found to be necessary for the 
dipole solitons at the interface of two nonlocal nonlinear 
media with a finite range of nonlocality. Here the ther- 
mal nonlinearity has an infinite range of nonlocality and 
the boundaries are essential [l8|. The boundary force 
effect, which vanishes in the nonlocal nonlinear media 
with a finite range of nonlocality, can support the sta- 
bility of dipole interface solitons, whatever small or large 
rid- From Fig. Etc), b c varies slightly when the sample 
size is doubled. It implies that the boundary force effect 
do not decrease significantly when the boundaries move 
far away from the solitons. Due to the infinite range 
of nonlocality, dipole interface solitons can exist in very 
large samples. 




FIG. 5. (Color online) (a) Propagation of a perturbed inter- 
face fundamental soliton at rid — 5, b — 0.5. Propagations of 
perturbed dipole interface solitons at (b) rid — 1, b = 16, (c) 
rid ~ 0.1, b — 1.5, and (d) rid = 0.2, b — 1.5. In all cases 
L = 30. 

To confirm the results of the linear stability analy- 
sis, we simulate the soliton propagations based on Eqs. 
(JTJ and © with the input condition q(x, z = 0) = 
ui(a;)[l + p(x)\, where w{x) is the profile of the station- 
ary wave and p{x) is a random function which stands for 
the input noise with the variance 5^ oise = 0.01. Figure[5] 
presents the propagations of fundamental and dipole in- 
terface solitons. As expected, the fundamental and dipole 
interface solitons at the stable region predicted by the lin- 
ear stability analysis survive over long propagation dis- 
tance in the presence of the input noise [Figs. Eta)-(c)]. 
Figure Efd) presents a dipole soliton in the instability 
region, it experiences oscillatory instability after propa- 
gating over a certain distance. 

To conclude, we have presented the study of interface 
solitons in thermal nonlinear media. This type of solitons 
shows some unique properties which differ from those of 
bulk solitons and surface solitons. The linear stability 
analysis shows that fundamental interface solitons are 
always stable and the stability of dipole interface soli- 
tons depends on the difference in the linear refractive 
index. We consider that the boundary force effect plays 
an important role in the stability of dipoles. As the index 
difference rid approaches to zero, our interface solitons re- 
duce to the corresponding bulk solitons which have been 
proved to be stable. For large rid, the interface solitons 
are similar to the surface solitons in intensity profiles, but 
different in nonlinear refractive indices. 
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